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ABSTRACT 


The  numerical  integration  technique  to  be  utilized  in  the  reduction  and 
analysis  of  data  gathered  in  the  USAF  Aeroballistic  Research  Facility  located 
at  Eglin  Air  Force  Base  is  described.  The  method  of  Chapman  and  Kirk  has 
been  developed  into  a  system  of  digital  computer  programs  which  may  be  easily 
and  routinely  used  in  the  analysis  of  data  derived  from  spark  range  firing. 

The  equations  of  motion  for  a  missile/projectile  in  free-flight  are  derived 
in  this  report  utilizing  the  six  degrees  of  freedom.  The  parametric  equations, 
philosophy,  and  methods  of  implementation  are  also  derived  and  discussed. 

Brief  descriptions  of  the  computer  program  involved  are  also  given. 
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SECTION  I 

INTRODUCTION  AND  BACKGROUND 


This  report  will  discuss  the  rationale  and  mathematics  involved  in  the 
implementation  of  a  numerical  integration  data  reduction  technique  to  be 
utilized  with  the  USAF  Aeroballistic  Research  Facility'1,  Air  Force  Armament 
Laboratory,  Eglin  Air  Force  Base,  Florida. 

2 

This  numerical  technique  was  reported  on  by  Chapman  and  Kirk  in  1969. 

3  4 

Goodman  (1966)  and  Knadler  (1969)  also  reported  on  this  general  technique 

*>6  7 

under  separate  efforts.  Whyte  and  Beliveau  of  General  Electric  Company  ’  * 

(1969,  1970)  applied  this  technique  to  ballistic  range,  yaw  sonde,  and  wind 

tunnel  data.  Beginning  in  1971  General  electric  Company,  under  various  in- 

house  efforts  and  United  States  Government  contracts,  has  developed  and 

applied  this  technique  to  a  wide  variety  of  exterior  ballistic  and  interior 

8-15 

ballistic  analysis  problems. 

In  October  1972  General  Electric  Company  delivered  to  Arnold  Engineering 

Development  Center  a  series  of  numerical  integration  data  reduction  computer 
13 

programs  for  use  with  VKF  Range  G.  These  programs  are  now  being  utilized 
in  the  routine  production  of  free-f light  data  from  Range  G. 

The  capability  exists  in  these  pregrams  to  analyze  up  to  three  experi¬ 
ments  simultaneously  and  to  identify  Mach  number  non-linearities  in  addition 
to  angle  of  attack  non-linearities. 

The  original  programs  were  obtained  by  AFATL  from  AEDC  and  have  been 
modified  to  complement  the  large  computational  capability  of  ADTC.  Most 
importantly,  the  programs  have  been  tailored  to  provide  maximum  flexibility 
of  use  while  decreasing  the  tasks  required  of  the  project  engineer. 

The  inputs  and  outputs  of  the  various  computer  programs  will  be  dis¬ 
cussed  briefly  in  this  report.  A  forthcoming  document  will  report  on  the 
details,  flow  charts,  input,  output,  and  options. 
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SECTION  II 

AEROBALLISTIC  RESEARCH  FACILITY 


The  United  States  Air  Force  Aeroballistic  Research  Facility  is  part  of 
the  Air  Force  Armament  Laboratory  located  at  Eglin  Air  Force  Base,  Florida. 
This  facility,  now  under  construction,  is  described  in  detail  in  reference  1. 
The  construction  is  expected  to  be  complete  and  instrumentation  installed  by 
summer  1974. 

This  facility,  essentially  a  ballistic  range,  is  constructed  of  concrete 
and  is  approximately  780  feet  in  length.  The  basic  instrumentation  will  be 
50  spark  shadowgraph  stations.  Each  station  will  consist  of  two  sparks  and 
two  cameras  located  in  the  wall  and  floor  (pit). 

The  stations  may  be  set  up  at  131  locations  along  the  range  thus  allow¬ 
ing  considerable  spacing  flexibility  to  enhance  data  acquisition  for  contem¬ 
plated  test  programs. 

Timing  data  will  be  provided  at  each  of  the  50  stations  by  electronic 
chronographs.  They  are  operated  with  an  IR  detection  system  and.  the  spark 
apparatus . 

The  estimated  accuracies  of  measurements  are  as  follows: 

Time  ±0.1  microsecond 
Position  ±0.0013  foot 
Angle  ±0.10  degree 

This  range  will  be  continuously  calibrated  through  the  use  of  catenary 
wires  suspended  in  front  of  the  wall  and  ceiling  reflective  screens.  Refer¬ 
ence  beads  will  be  precisely  positioned  along  the  wires,  with  the  location  of 
all  beads  known  to  0.0001  foot. 
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SECTION  III 
STATEMENT  OF  PROBLEM 

The  angular  and  translational  motion  of  a  projectile  as  it  traverses  a 
ballistic  range  may  be  observed,  measured,  and  documented  by  means  discussed 
in  the  previous  section  and  reference  1.  The  following  basic  data  are  ob¬ 
tained  as  a  function  of  time:  X  ,  Y  ,  Z  ,  0,  tj;,  and  <j>.  Knowing  the  physical 
properties  of  the  projectile  and  the  atmospheric  properties  of  the  range 
facility,  the  problem  becomes  one  of  identifying  the  aerodynamic  coefficients 
and  their  magnitudes  which  are  causing  the  observed  motion. 

The  most  prevalent  method  of  analyzing  ballistic  range  data  is  known  as 

linear  theory.  Murphy***  MacAllister*^’^,  and  others  of  BRL  as  well  as 
19  21  22  23 

Nicolaides  1  ’  and  Eikenberry  have  developed  this  theory  to  the  extent 

that  certain  types  of  non-linearities  car  also  be  analyzed.  However,  unless 
many  cycles  of  data  are  present  and/or  multiple  experiments  conducted,  the 
only  non-linear  coefficient  which  can  be  consistently  identified  is  the  cubic 
pitching  moment. 

In  1969  Chapman  and  Kirk  of  NASA  Ames,  in  attempting  to  analyze  free 
flight  data,  developed  a  technique  which  allowed  the  differential  equations 
of  motion  to  be  used  directly  in  the  fitting  process. 

Linear  theory,  on  the  other  hand,  depends  on  using  the  closed  form 
solution  to  the  equations  of  motion  in  the  fitting  process.  This  technique 
is  very  effective  providing  the  assumptions  required  to  obtain  the  closed 
form  solution,  i.e.,  small  velocity  drop,  nearly  constant  pd/2V,  linear  aero¬ 
dynamics,  and  small  angular  motion  are  not  drastically  violated. 

None  of  the  above  assumptions  are  required  when  using  the  differential 
equations  of  motion  directly.  The  only  drawbacks  or  restrictions  are  in  the 
identification  of  the  coefficients  which  are  to  be  solved  for  and  the  efficient 
coding  of  a  computer  program  such  that  a  solution  may  be  found  at  a  reasonable 
cost. 
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The  system  in  use  at  AEDC  and  Eglin  AFB  uses  the  linear  theory  programs 
in  concert  with  the  numerical  integration  programs  such  that  .the  user  has 
fewer  problems  in  obtaining  adequate  reductions.  Linear  theory  reduction  and 
numerical  integration  reductions  may  be  directly  compared  as  they  occur  in  the 
same  computer  run. 
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SECTION  IV 
THEORY  AND  EQUATIONS 


1.  EQUATIONS  OF  MOTION 

The  equations  of  motion  are  derived  in  a  fixed  plane  coordinate  system. 

The  Euler  .ingles,  as  well  as  the  missile-fixed  velocities  and  earth-fixed 
velocities,  are  integrated  numerically  (step-wise)  to  simulate  the  projectile 
motion. 

This  formulation  was  chosen  for  the  following  reasons: 

1.  Allows  largest  time  step  for  integration. 

2.  Most  accurate  for  spin-stabilized  or  fast  spinning  projectiles. 

3.  Consistent  with  data  acquisition. 

4.  Allows  slight  aerodynamic  asymmetries. 

The  following  assumptions  are  contained  in  the  derivations: 

1.  No  earth  rotation  effect  on  angular  motion 

2.  Flat  earth 

3.  Rotational  symmetry  (C  =  f.  I  =  I  ,  etc.) 

met  Up  y  z 

4.  Trim  moments  and  forces  are  small. 

5.  Aerodynamic  coefficients  are  expanded  as  polynomial  functions 
of  the  sine  of  the  total  angle  of  attack. 

These  derivations  are  patterned  after  Barnett's  of  Picatinny  Arsenal.2^’25 
These  references  are  interesting  in  that  they  describe  the  advantages  of  one 
type  of  coordinate  system  as  opposed  to  another,  Barnett  discusses  the  pro¬ 
blem  associated  with  obtaining  adequate  (accurate)  simulations  of  spin  sta¬ 
bilized  projectiles  using  a  coordinate  system  which  is  attached  to  the  projec¬ 
tile  and  rolling  with  it.  His  conclusions  and  those  of  the  authors  of  this 
report  are  that  adequate  simulations  of  spin  stabilized  projectiles  may  most 
reliably  be  made  with  the  fixed  plane  system.  In  addition  a  run  time  cost 
savings  of  over  25  to  1  are  realized  when  using  a  fixed  plane  simulation  as 
opposed  to  a  rolling  coordinate  system  simulation  of  questionable  accuracy. 
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a.  Six  Degrees  of  Freedom 


Force  Equations  of  Motion: 


(Earth-fixed  derivative) 


CD 


l  F  =  Summation  of  forces  acting  on  missile 

R  =  Vector  from  earth  center  to  current  missile  C.G. 

i  ,  j  ,  F„,  =  Unit  vectors  aligned  with  missile- fixed  coordinate  system 
H  H  (non-rotating) 

a)  =  Angular  velocity  of  missile- fixed  coordinate  system  relative  to 
earth- fixed  coordinate. 


6 


(5) 


where:  V  =  u  i„  +  v  i„  +  w  k„ 
n  ri  rl 


dHV  -r 

d  t 


(6) 


where  all  quantities  in  the  right-hand  side  of  Equation  (4)  are  understood 
to  be  expressed  in  missile  coordinates  (H) .  Intuitively,  the  presence  of  the 
a)  x  V  term  relates  the  motion  of  the  missile  coordinate  system,  to  which 
forces  and  moments  are  referred,  to  earth-fixed  coordinates. 

Combining  Equations  (4) ,  (5) ,  and  (6) ,  and  performing  the  indicated 
operations,  one  then  obtains  the  three  component  force  equations: 


E  FXH  =  m  [u  +  (uYH  w  -  o:ZH  v};(in  the  i  direction} 
♦ 

£  FYH  *  ■  Cv  *  “ZH  U  -  “XH  W 

z  fzh  =  m  0  *  “XH  V  -  “YH  u^»{in  the  direction} 


Moment  Equations  of  Motion: 

Treating  the  moment  equation  in  a  similar  manner  to  the  force  equation: 

—  dET  “h"7 _ _ 

l  L  =  -  =  -  +  a)  x  J  (8) 

d  t  d  t 

Summation  of  moments  acting  on  missile 
Angular  momentum  of  missile 

Total  angular  velocity  of  the  missile  (rolling  coordinates) 
relative  to  earth  coordinates. 


}»(in  the  jH  direction} 


(7) 


(9) 


(oj^xh  XH  +  \  (“■pYH  +  ly  ^VzH  kH 


where  I  ,  I  are  the  moments  of  inertia  of  the  missile  about  the  longitudinal 
x  y 

and  traverse  axes,  respectively.  Note  that  the  assumption  of  rotational 
symmetry  (I  =  I  )  is  implicitly  made. 

z  y 


v  .  „  .  . 

"  ■  *x  ^VxH  XH  *  :y  ^VyH  *  Ty 

d  t 


ZH 


'H 


CIO) 


Appropriately  combining  Equations  (2) ,  (.8) ,  (9) ,  and  (10)  results  in 
the  component  equations  for  the  moments: 


E  LXH  =  !x  (VxH  +  WYH 
1  =  Iy  ^VyH  +  “ZH 
1  LZH  =  Ty  (<VzH  +-C0XH 


ly  ^VzH  "  WZH  ly 
*x  ^VxH  "  WXH  ly  (VzH 
ry  (“T-’YH  "  “YH  Ix  ‘■VxH 


(11) 


AERODYNAMIC  FORCES 


Three  aerodynamic  forces  are  considered: 

Axial  drag  -  Along  missile  axis 

Normal  force  -  Perpendicular  to  missile  axis  in  plane  of  yaw 

Magnus  force  -  Perpendicular  to  plane  of  yaw 


These  forces  are  illustrated  in  Figure  1  (arrows  indicate  positive 
directions) . 


Figure  1.  Force  System 

The  axial  drag  acts  along  the  negative  axis  which,  by  definition,  is 
directed  from  the  CG  towards  the  nose  of  the  missile. 


(Axial  Drag) m  = 


p  d2  it  (V)2 


-  q  A  C, 


Where  q  =  dynamic  pressure  =  1/2  p  v 


A  =  cross  sectional  area 


The  normal  force  components  act  opposite  to  the  directions  of  v  and  w. 
Using  known  geometry,  one  can  deduce  the  components  of  the  normal  force: 


(Normal  Force) Y„  =  -  q  A  Cw  sin  a 


'YH  -  ~  "  "No  "***  “  J2  2 

Vv  +  w 

where:  =  The  slope  of  the  Normal  Force  coefficient  plotted  as  a  function 

of  sin  a. 
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.*  (Normal  Force)yH  =  -  q  A  CNa  I 


(Normal  Force) ^  =  -  q  A 


sin  a 


w 

V 


The  Magnus  Force  acts  in  a  direction  perpendicular  to  the  plane  of  yaw. 
As  in  the  case  of  the  Normal  Force,  using  known  velocity  geometry  the  Magnus 
Force  components  are: 


(Magnus  Force)  m  (VxH  CYpa  sin  ° 


2  2 
v  +  w 


Where  Cypa  =  The  slope  of  the  Magnus  Force  coefficient  plotted  as  a  function 
of  sin  a. 

(Magnus  Force) w  =  S_*-i  (^  S 
(Magnus  Force)  ZH  =  -  q2  v  (VxH  ®Ypo  V 


Summary  of  Forces:  (excluding  gravity) 

1  FXH  =  "  ^  A  S 

1  FYH  =  "  **  A  Sta  v/V  +  ^2  V  XII  Sfpa  V  ^ 

1  FZH  =  “  *  A  Sta  V  "  ^2  V  '  ^XH  Spa  V 
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AERODYNAMIC  MOMENTS 


The  following  moments  are  considered: 


Pitching  Moment 
Magnus  Moment 
Damping  Moment 
Spin  Deceleration 
Trim  Moment 


-  Perpendicular  to  plane  of  yaw 

-  In  plane  of  yaw 

-  Perpendicular  to  missile  axis 

-  Along  missile  axis 

-  Perpendicular  to  missile  axis 


Rotational  sign  convention  is  indicated  in  Figure  1. 


Pitching  Moment  Components: 


(Pitching  Moment) v„  =  q  A  d  C  sin  a 
in  ma 


Where:  *  The  slope  of  the  Pitching  Moment  coefficient  plotted  as  function 

of  sin  a. 


(Pitching  Moment) YH  =  q  A  d  ~ 

(Pitching  Moment)  =  -  q  AdC 

4ri  ma  V 


Magnus  Moment  Components: 


(Magnus  Moment) YH 


q  A  d2 
2  V 


CVxH  Cnpa  Sin  a 


Where:  C  o  =  The  slope  of  the  Magnus  Moment  coefficient  plotted  as  a  function 
of  sin  a*. 


(Magnus  Moment)^  .  c„p„  V 

(Magnus  Moment) Cnpa  H. 
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Damping  Moment  Components: 


-  A  d2 

(Damping  Moment) ^  —  (wT)yH  C 


mq 


Where:  C  =  Damning  Moment  coefficient, 
mq 


“Ad2 

(Damping  Moment)  m  =  —  (o>T)ZH  C 


mq 


Spin  Deceleration  Moment: 


“Ad2 

(Spin  Moment)^  =  3__  (^  CJp 


=  Spin  deceleration  coefficient 


Trim  Moment  Components: 


This  moment  is  induced  by  missile  asymmetries  and  is  assumed 
for  spin- stabilized  missiles. 


(Trim  Moment)^  =  -  q  A  d  [C^  (6g  sin  <f>  -  6^  cos  4>)  ] 


(Trim  Moment)  =  q  A  d  [Cm<j  (6g  cos  <f>  +  6^  sin  <j>)  ] 

Where:  C  .  =  Trim  Moment  Coefficient 

mo 

6.,  =  Components  of  the  trim  angle 

A  D 

<|>  =  Roll  angle 


negligible 
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Moments : 


1  a  i_ 

2  V 


‘•VxH  Clp 


I 

x 


^VxH  +  “YH  ly  ^VzH 


ID-u  I 
ZH  y 


q 


A  d  C 


ma 


q  A  d2 
2  V 


O^XH 


C 

npa 


v  + 
V 


q  A  d2 
2  V 


^Vyh 


C 

mq 


q  A  d  [Cm5  C6B  Sin  *  -  6A  005  *)]  =  \  CVYH  +  WZH  Tx  (VxH 


“XH  ly  (VZH  1  (nation  jH  j 


(18) 


—  2  —  2 
—  .  ,  7T  v  q  A  d  ,  .  pr  w  q  A  d  .  .  77 

q  A  d  Snct  V  H  UT  XH  Cnpa  V  2  V  ^ZH  Snq 


+  q  A  d  [Cm6  (*B  C0S  ♦  +  6A  Sin  "  Xy  (VZH  +  “XH  Xy  CVyH 


UYH  Xx  (VxH  ; 


rotation  k, 


i.) 


(19) 


Equations  (14)  through  (19)  are  of  limited  use  until  one  knows  explicitly 
how  the  missile -fixed  coordinate  system  (H)  is  moving  with  respect  to  earth- 
fixed  coordinates  (E) .  Further,  one  often  has  vector  quantities  expressed  in 
earth  coordinates  (such  as  wind  and  gravitational  attraction)  which  must  be 
properly  introduced  into  the  (H)  coordinates.  Therefore,  additional  relation¬ 
ships  between  the  (H)  and  (E)  coordinates  must  be  derived.  In  particular,  one 
must  know  how  the  (H)  coordinate  system  is  oriented  relative  to  the  (E)  coordi¬ 
nates  at  all  times.  The  use  of  Euler  angles  seems  to  be  the  most  straight¬ 
forward  approach.  In  it,  one  rotates  a  coordinate  system,  initially  coinci¬ 
dent  with  the  (E)  system,  about  selected  axes  so  that  after  the  rotations 
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are  performed  in  a  specified  sequence,  this  coordinate  system  will  have  the 
same  orientation  as  the  (H)  system.  These  rotations  are  illustrated  by 
Figure  2. 


Z 


x,  y,  z  -  Earth  Axis  Coordinate! 

x',  y',  z'  -  Intermediate  Axis  Coordinates 
(after  first  rotation  -  through  'F) 

X,  Y,  Z  -  Missile  Axis  Coordinates 

6  and  'I'  are  the  Euler  angles 


Figure  2.  Euler  Coordinate  System 


Figure  2  shows  the  (H)  and  (E)  coordinate  system.  The  missile-fixed 
rolling  coordinate  system  (M)  is  defined  as  simply  the  (H)  coordinate  system 
rotating  about  the  missile  axis  at  an  angular  rate  of  <J>  through  the  angle  <j>. 
Previously  defined  quantities  in  the  equations  of  motion  may  now  be  related 
through  the  Euler  angles  to  earth-fixed  coordinates  as  follows: 


^VxH 

cVyh 

(VzH 

(VXH 


Wjjjj  =  -  sin  0 

“YH  =  6 
wZH  =  \p  cos  e 


UXH  =  -  “ZH  tan  6 


$  +  =  <}>  -  ifi  sin  e 


(20) 


UYH  ° 


a>ZH  =  ip  cos  0 


<f>  -  uZH  tan  0  -  u>ZH  6  sec  0 

•  •  .♦  •  •  •  • 

<J>  -  tan  0  (ij»  cos  0-^/0  sin  0)  -  0  ij)/cos  0 


•  •• 

(w^yH  =  6 

•  ••  •  • 

(o)T)  ZH  =  cos  0  -  ip  0  sin  0 


Through  appropriate  substitution  of  Equations  (20) ,  the  equations  of 
motion  may  be  rewritten  as : 

Force  Equations: 


.-qACy  =  m  [u,+  Ok  -  ^  cos  0  v] 


;  jjn  iH 


(21) 

(X)  direction^ 
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-  -  v  ,  q  A  a  w  «i  .  4,  sin  0)  *  »'  W  *  *  “  cos  0 

q  A  Ckl_  „  +  9  v  '‘Ynix  V 

S In  JH  (Y)  direction 


-  H  A  '■‘Net  V  2  V  Yp«  V 

+  w  sin  0] 


_  .  -  w  q  A  d  r  I  (i  -  i  sin  0)  =  m  [w  -  i  v  sin  6  (23) 

-  q  A  cNa  v  -  2  V  Ypa  V 


-  0  u] 


j  fin  kn  (Z)  direction 


Moment  Equations: 


(i  .  ;  si,  e)  -  I.  li  -  tan  e  (i  cos  9  -  i  e  sin  e)  (24) 

2V  l’lp  X 


-  iji  0/cos  0] 


;  ^Rotation  about  (X)  ^ 


-  •  .2  «  ‘  •  -  +  q  A  d"  q  c  _  (25) 


-  A  ,  7f  H  +  S  A  j-  C  J-  C4>  -  *  sin  0)  +  Vv  '  _mq 

*1  A  d  Cmo  V  2  V  *Pa  V 


-  q  A  d  [Cm6  (6b  sin  $  -  cos  *  Xy 


•  » 

X  e  +  I  4»  cos  0  (<t»  -  sin0) 


.2 


+  I  it  cos  0  sin  0 

y 


:  /  Rotation  about 

l  * — 


-  d  5  V  add*  *  X.  (i  -  *  sin  «  *  Vt-  cmq  * cos  9  (26) 
•  ‘i  A  “  lia  V  2  V  npa  V 


♦  ?  A  d  [C.j  C«B  “S  d  +  SA  si”  ’  V  C*  C°S  9  ’  *  6  S1"  6) 


•  *  ,  *  t 1  i  *  (  Rotation  about  k„  (Z) 

♦  !  C-  *  *in  8)  e  -  Ix  9  (♦  “  *  sm  0)  ’  '  H 


Equations  (21)  through  (26)  are  the  6-DOF  equations  of  motion  expressed 
in  terms  of  Euler  angles  without  gravitational  effects. 

Gravitational  Force: 

Assumption:  Flat  Earth 
(Gravity)  =  -mg  Fg 

where  F_  is  unit  vector  in  the  +ZC  direction 
E  E 

Replacing  Fg  by  its  representation  in  (H)  coordinates,  one  obtains: 

-mg  kg  =  mg  sin  6  TH  -  mg  cos  6  kj^ 

The  gravity  force  can  now  be  substituted  directly  into  Equations  (21) 
and  (23) . 

Spin  is  defined  as  the  total  angular  velocity  of  the  missile  about  its 
own  axis. 

•  • 

Spin  =  p  =  (u^)^  sin  8 


The  6-DOF  equations  of  motion,  upon  substitution  of  the  gravitational 
force  and  spin  identity  and  rearranging  terms,  are  written  as  follows: 


Force  Equations: 


u 


«  S  cx  *  * 


•  » 

sin  8  -  0  w  +  ij>  cos  0  v 


(27) 


-  A  r  v  . 
'  q  m  Net  V 


q  A  d  p  -  w 
m  2  V  Ypa  V 


-  cos  0  [u  +  w  tan  8] 


(28) 
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(29) 


w  =.  -  a J-r  £- 

m  Na  V 


q  A  pd  p  v 

m  2  V  Ypa  V  " 


g  cos  0  +  i{)  (sin  0)  v  4  0  u 


Moment  Equations: 


-  .  j2 


••  Q  A  d  D  **  •  *  •  • 

*  *  I„  2V  P  Clp  +  tan  0 tip  cos  0  -  0  sin  0]  4  e  ij;/cos  0 


(30) 


0  =  q. A  d  c  -  4  q  A  d2  c  0  4  q  A  d2  P  r  v 
ly  ma  V  2V  I  mq  8  +  2V  I  Sipa  V 


(31) 


[Cm6^6B  Sin  *  "  6 A  cos 


i]>  cos  9  p  j— 


.2 

cos  0  sin  0 


-  .  j2 


i  -  [  _  q  A  d  C  I  4  q  A  d  c  i  rn,  fi  4  q  A  d  p  p  W 

*  I  ma  V  2V  I  Cmq  *  C0S  6  Yv"  Iv  Cnpa  V 


.(32) 


T  i  j  •  # 

+  iL -  [Cmg  (6g  cos  4>  4  6  A  sin  <J>)]  4  2  0  i|>  sin  0 


4  0  p  I  /I 

r  x  y 


]  /cos  0 


For  the  data  reduction  technique  being  discussed,  the  aerodynamic  coeffi¬ 
cients  C  ,  C  ,  and  C  previously  defined  are  assumed  to  be  non-linear  and 
ma'  npa'  mq  r 

are  therefore  expanded  to  include  higher  order  dependencies  on  the  angle  of 
attack,  a.  These  dependencies  are  assumed  to  exist  as  a  function  of  sin  a. 
The  following  non-linearities  are  assumed: 
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—  2  —  4  —  6  — 

C  =  C  +  C  sin  a  +  C  sin  a  +  C  sin  a  +  C 
ma  ma  ma3  mas  mct7  mot 


<Vo 


-  V) 


C  =  C  +  G 
mq  mq  mq?_ 


.  2  - 
sm  a 


C  =  c  +  c 
npa  npa  npa  3 


sin"  a  + 


C 

npa  5 


.  4  - 
sm  a 


The  preceding  moment  coefficient  expansions  are  utilized  in  the  angular 
motion  programs. 

Similar  expansions  on  the  force  coefficients  are  used  in  the  translational 
equations  of  motion  which  will  be  developed  next. 

b.  Translational  Equations  of  Motion 

Tn  the  derivation,  the  previously  defined  force  equations  of  motion  will 
be  utilized  with  the  addition  of  a  trim  force  in  the  missile  Y-Z  plane,  de¬ 
fined  as  components  in  the  missile  Y  and  Z  directions.  This  trim  force  will 
account  for  any  slight  asymmetries  present  (intentional  or  otherwise)  in  the 

missile  configuration. 

•  »  •  • 

u  =  -  3—  Cv  +  g  sin  0  -  0  w  +  <|/  cos  0  v  (33) 

m  a 


=  .  ±J±C  _  +  ii_£l  1L2±  c 

m  Na  V  m  2  V  Ypa 


0  [u  +  w  tan  0] 


-  V  tS8  C0S  *  '  5A  Sin  ♦> 

A 


(34) 


W  =  -  SLA  c  _  lA  C  ^  -  g  cos  0  +  |(|  sin  0  v 
w  m  Na  V  m2V  Ypa  V 


+  0  u  -  CN($  (<$A  cos  <j>  +  <5g  sin  <{>) 
A 


/  _  ) 
q.A 
m 

j 


(35) 


20 


The  following  expressions  relate  the  earth-fixed  velocities  to  the  missile 
velocities  using  the  definitions  of  the  Euler  angles: 


V^e  =  u  cos  0  cos  if  -  v  sin  if  +  w  sin  G  cos  if 
Vy  =  u  cos  0  sin  if  +  v  cos  if  +  w  sin  0  sin  if 


=  -  u  sin  0  +  w  cos  0 


(36) 


Taking  the  derivative  of  Equations  (36)  With  respect  to  time,  one  obtains: 
•  *  *  • 

VXe  =  u  cos  0  cos  if  -  v  sin  if  +  w  sin  0  cos  if  (37) 

♦  •  • 

-  u  sin  0  cos  if  0  -  u  cos  0  sin  if  if  -  v  cos  if  if 

*  • 

+  w  cos  0  cos  if  0  -  w  sin  0  sin  if  if 

•  •  »  « 

Vye  =  u  cos  0  sin  if  +  v  cos  if  +  w  sin  0  sin  if  (38) 

♦  •  • 

-  u  sin  0  sin  if  0  +  u  cos  0  cos  if  if  -  v  sin  if  if 

•  • 

+  w  cos  0  sin  if  0  +  w  sin  0  cos  if  if 


VZe  =  -  u  sin  0  +  w  cos  0  -  u  cos  0  0  -  w  sin  0  0  (39) 


Substitution  of  Equations  (33) ,  (34) ,  and  (35)  into  (37) ,  (38)  ,  and  (39) 
provides  the  translational  equations  of  motion: 


plAv2 

2  m 


cos 


0  COS  If 


*f£v2 


CNa  V  si"  * 


(40) 


.  lA  v2 

2  m 


sin  0  cos  if 


p  A  pd 
4  m 


7T  W 


S 


pa 


sin  if 
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p  A  pd 
4  m 


—  v 

V  C„  rr  sm  0  COS 
Ypa  V 


P 


+ 


p_A  v2 

2  m 


[CNg  (<Sg  cos  <f>  -  $A  sin  4>)  sin  p 
A 


"N6a  (6  cos 


4  +  <$, 


sin  <f>)  sin  8  cos  ifi] 


.lAv2 

'Ye  2  m  "X 


V.,  =  -  Cv  cos  6  sin  \J»  -  —A  V2  C..  cos  p 

Y  2  m  Not  V 


(41) 


p  A  ,,2  rr  w  .  A  .  ,  ,  p  A  pd  ,,  77  w 

£ -  V  CL,  rr  sm  0  sm  ^  V  C„  rr  COS  if/ 

2  m  wx  V  4  m  Ypa  V 


p/-  V  Cv  ^  sin  8  sin  p 
4  m  Ypa  V  r 


p  2 

2*^  V  [CN6  (6g  cos  p  -  6a  sin  <j>)  cos  i|» 
A 


+  CN5  (6a  cos  4>  +  6g  sin  <j>)  sin  0  sin  i|j] 
A 


V7  =  2-A  V2  CY  sin  8  -  V2  CM  £  cos  0 
Ze  2  m  X  2  m  Na  V 


(42) 


-  V  Cv  w  cos  0 

4  m  Ypa  V 


'n6  (<5a  cos  4»  +  6g  sin  <j>)  cos  8 

A 


£.  v2 

m 


The  following  non-linearities  are  assumed  for  the  force  coefficients: 

ex  =  cx  *  %  \  (V  -  V0) 

2 


C.,  =  C,  +  CM  e 
Na  Na  Na3 
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Ypa  ipa  Ypas 


2  2  - 

Where:  e  =  sin  a 

VQ  =  Reference  Velocity  (constant) 


c.  Modified  Six  Degrees  of  Freedom 

As  ballistic  ranges  are  level  and  the  ballistic  trajectories  flat,  the 
equations  of  motion  from  which  the  parametric  equations  are  derived  can  be 
modified.  This  is  a  special  case  and  cannot  be  applied  to  high  angle  fire 
problems  (y  >  30°)  such  as  yaw  sonde  reductions.  Reference  11  discusses 
the  consequences  of  this  assumption  in  the  high  angle  fire  mode. 

The  following  assumptions  are  made: 


and  from  linear  theory  (Ref.  17) 


C  =  C  -  C.r  2  T  /md 

mq(3D)  m<i  N«  y 


:  =  C  +C„  2  1  /md 

nPa(3D)  npa  Na  x 


Essentially  what  is  done  is  to  approximate  the  6-DOF  equations  with  3-DOF 
equations  and  allow  correction  for  the  effect  on  the  damping  of  C^. 

The  effect  of  trajectory  curvature  is  accounted  for  by  the  angular  rate 
g  cos  e/v. 


9 (3D)  =  °(6D) 


+  g  cos  e/v 


23 


This  term  results  in  the  generation  of  the  yaw  of  repose. 

It  should  be  pointed  out  and  understood  that  the  modified  equations  are 
only  used  in  deriving  the  parametric  equations  and  are  not  used  for  computing 
theoretical  motions.  Equations  (28),  (29),  (31),  and  (32)  are  utilized  for 
motion  generation. 

Approximate  6-DOF  Equations 

Using  the  previous  assumptions,  Equations  (31)  and  (32)  reduce  to: 

••  * 

6  =  M  sin  0  cos  ip  +  M  [0  +  g  cos  6/V] 
ma  mq 


M  sin  ip  -  M  . 
npa  Y  m6 

2 


0 


f  cos  0  p  Ix/I 


-  ip  cos  0  sin  0 


t 


t|)  «  [M  sin  <J>  +  M  ip  cos  0  +  M  sin  0  cos  ^  +  M. 


ma 


mq 


npa 


m6 


I  .  . 

+  (0  +  g  cos  0/V)  p  +  2  sin  0  (0  +  g  cos  0/V)  i^]/cos  0 

y 


Where: 


M  =q~[C  +  C  c 

ma  n  I  ma  ma.-? 


2  +  C  e'  +  C  (V„  -  V)] 
ma3  mas 


may  o 


-  .  .2 


M  =  3-44-  [C  +  C  e  -  CM  2  I  /md2] 
mq  2  V  I  mq  mq2  Na  y 


-  .  ,2 


M  =  3  rr  T--  p[c  +c  G+C  e+Cw2  I  /md2] 
npa  2  V  I..  v  1  npa  npa,  npas  Na  x 


npa  npa3  npas 
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MmS  *  1  f1  tCm6  ({B  sin  *  '  SA  cos  W1 
e  y 


V  *  ’  r1  [C«S  tSB  cos  +  *  5A  si"  «' 

♦  y 


As  the  term  g  cos  0/V  is  usually  small  and  1.0  >  cos  0  >  0.9,  this  term 
is  considered  a  constant  when  partial  derivatives  are  taken. 

d.  Roll  Equation 

As  previously  derived  and  stated  in  Equation  (30)  the  6-DOF  roll  equation 
is: 


; . 11 A 


VIX  -ip 


C1  (4>  -  if)  sin  0)  +  tan  0  (if)  cos  0  -  if;  0  sin  0)  + 


if;  0 
cos  6 


The  preceding  equation  may  be  restated  in  the  following  equivalent  form: 


#  =  p  +  if)  sin  0 


(30a) 


•  q  A  d2 
p  -  %^-j-  p  C 


IP 


(30b) 


The  above  two  equations  are  integrated  simultaneously  to  produce  the  roll 
angle,  roll  rate,,  and  spin  rate  with  the  advantage  of  fewer  higher  order  terms 
than  equation  (30) . 


if)  will  be  computed  from  the  available  if)  versus  time  history.  It  should 
be  pointed  out  that  the  roll  reduction  is  accomplished  at  the  same  time  as 
the  translational  motion  reduction. 

—Ad  r 

A  fin  cant  moment  (Cj5  6  q  -j — •)  and  a  velocity  dependent  C 
—  i  ,2  x  <— 


<V0  -  V>  2  V 


—  2 
SL  ft  il£ 

2  VIx  J 


lpv 

are  added  to  equation  (30b)  for  expanded  capability. 
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2.  CHAPMAN-KIRK  TECHNIQUE 


The  procedures  utilized  in  the  implementation  of  the  Chapman-Kirk  technique 
are  as  follows: 

1.  Formulate  equations  of  motion  which  are  adequate  to  simulate  the  ex¬ 
perimental  data. 

a.  Select  a  broad  general  set  of  aerodynamic  terms. 

b.  Define  possible  non-linear  terms;  i.e.,  Mach  number  and  angle  of 

attack. 

2.  Equations  of  motion  are  partially  differentiated  with  respect  to  each 
coefficient  to  form  a  set  of  parametric  differential  equations. 

3.  Integrate  equations  of  motion  numerically  utilizing: 

a.  Estimated  Aerodynamic  Coefficients 

b.  Estimated  Initial  Conditions 

4.  Integrate  parametric  equations  numerically  to  obtain  values  for  the 
partial  derivatives  of  the  state  variables  with  respect  to  each  coefficient. 

5.  A  differential  corrections  equation  is  set  up  from  a  Taylor  expansion 
of  the  dependent  variable  of  the  equations  of  motion. 

6.  The  experimental  data  are  then  compared  to  the  dependent  variables  in 
a  least  squares  sense»and  corrections  to  the  coefficients  and  initial  condi¬ 
tions  are  obtained  from  the  differential  correction  equations. 

7.  Step  3  is  repeated  with  adjusted  coefficients  and  initial  conditions, 
and  the  process  is  repeated  until  convergence  is  achieved. 


a.  Differential  Corrections  and  Least  Squares  Theory 


The  methods  of  least  squares  and  differential  corrections  are  employed 
by  computer  programs  to  obtain  corrections  to  be  applied  to  the  coefficients 
so  that  the  solutions  to  the  equations  of  motion  are  a  better  fit  to  the  test 
data.  A  brief  description  of  this  technique  will  be  given  considering  two 
data  sets.  The  primary  goal  of  this  description  is  to  illustrate  how  two  or 
more  data  sets  may  be  handled  simultaneously.  Consider  the  following 
equations: 

al  +  ^lal  +  <"'2al  +  ^3  =  0 
•  •  • 

a2  +  Cla2  +  °2a2  +  C3  =  0 


=  Angles  from  data  set  one 


a2  =  Angles  from  data  set  two 

Although  the  above  equations  are  linear,  it  is  well  known  that  aerodynamic 
equations  of  motion  are  often  non-linear  in  nature.  Regardless  of  the  linear¬ 
ity  or  non-linearity  of  the  equations  of  motion,  the  theory  of  least  squares 
is  not  sufficient  in  itself  to  solve  for  aerodynamic  properties,  and  differ¬ 
ential  corrections  are  employed.  The  above  equations  are  presented  only  for 
simplicity  and  clarity.  The  real  equations  of  motion  cover  all  six  degrees 
of  freedom  and  are  coupled. 


Let: 


=  initial  condition  on  ct^ 


•  « 
Cj.  =  cijQ  =  initial  condition  on 


C 


6 


o2q  =  initial  condition  on  a2 


C?  =  a2g  =  initial  condition  on  a2 


Cj,  C2,  and  are  constant  coefficients  to  be  determined  and  are  common 
to  both  data  sets. 
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Let: 


t  =  The  variable  parameter  (time) 


<*1  =  f  CC,,  &2>  ^3»  ^4»  *) 

^2  =  8  (Cj  >  C2  y  >  Cg  t  Cy ,  t) 

Eij  =  f  (C1»  C2*  C3*  C4*  C5*  V  "  alj 
E2 j  =  g  (C1’  C2  *  C3  *  C6’  C7’  tj)  "  a2j 


J 

•  “2' 

E_ . 

n 

E 

-  •  ai. 

ij  ’  *2i 

computed 

A 

L.  .  a.  . 

b. 

►  experimental 

- ^ 

t. 

J 


With  first  approximations  of  the  constants,  a  Taylor  series  expansion  gives: 


+  Ejj  =  f  (Cj ,  C2,  Cj,  C4,  Cg,  tj)  ~  f  (C^,  tj)  +  i  Cj  f^ 


+  A  C,  f  +  A  C_  f  +  AC,  f 
2  c2  3  c3  4  c4 

[AC,]2  [AC]2 

+  A  C-  f  +  XT’ —  f  +  — ST -  f 

5  Cr  2!  c,  2!  c9 

b  l2  l2 

[A  C-]2  [A  CJ2  [A  C,]2 

+  -ttl-£c3  +ir"fc.  +-2r-£c5  +  •  •  • 
32  42  52 


°2j  +  E2j  =  g  (C1*  C2’  C3’  C6*  CV  V  3  g  CCo'  V  +  h  C1  gc, 


+  A  C2  gc  +  A  C3  gc  +  A  C4  gc 
2  3  4 
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+  A  Cc  g-  + 
*  C5 


[A  Cj]' 


[A  c2r 


5  gc_  +  2!  gc,  +  2!  gc„ 


[A  C  ]2  [A  C412  [A  C5]2 

+  —2!  gc  +  —J?  gc  +  — 21  gc  +  ' 

32  62  7 2 


.+  (higher  order 
terms) 


Where: 


al j ’  a2j  are  exPerimenta^  P°ints  of  data  sets  one  and  two,  respectively. 

E. .  E„ .  are  the  errors  between  the  data  and  the  theoretical  points  of 
data  sets  one  and  twc*  respectively. 


f  (CQ,  t)  =  f  [initial  guesses  at  constants  and  initial  conditions  (cij) ,  t) 

g  (CQ,  t)  =  g  [initial  guesses  at  constants  and  initial  conditions  (a2) >  t] 

A  CL  is  the  change  in  (L  between  successive  approximations. 

are  the  partials  of  f  and  g,  respectively,  with  respect  to  C 

c.’  gc. 

1  x 


f  ,  g  are  the  second  partials  of  f  and  g,  respectively,  with  respect 
x2  2  to  CL 

Assuming  that  the  A  C's  are  small  such  that  the  second  and  higher  order 
terms  can  be  neglected,  the  following  set  of  equations  result: 

“lj  ♦  *1J  -  f  V  ♦  »  ct  fCj  *  4  c2  fCj  uc,fCj(1  C4  fC4  *  »  c5  fCs 
°2j  *  E2j  c  6  CC0,  tj)  C2  *  4  C;  (Cj  *  1  Cg  g^  *  4  C, 


Let  the  residuals  be  the  error  at  each  point  that  result  from  using  the 
initial  guesses  of  the  constants  and  initial  conditions. 
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E  E1  £cx  *  E2  *Cj 

£  E1  %  *  E2  eC2 

2  E1  +  E2  gc 

1  c3  z  c3 

1  E1  fcA  +  E2  gcA 

4  4 

Z  E1  £c.  +  E2  % 

J  J 

1  E1  fc,  *  E2  *<• 

6  6 

£  E1  fc7  +  E2  8c7 


[A  C] 


The  summation  is  from  1  to  N,  where  N  is  the  number  of  data  points  in 
each  data  set. 


f  ,  f  ,  g  ,  g  do  not  appear  in  the  equations  and  are  equal  to  zero. 
c6  c7  c4  c5 

The  reason  for  this  is  that  the  initial  conditions  are  unique  to  each  data 
set,  as  opposed  to  the  common  set  of  coefficients  used  by  all  data  sets. 

b .  Parametric  Equations 


(1)  Derivation  of  Generalized  Parametric  Differential  Equation  for 
Translational  Motion 


Consider  equation  (40) : 


'xe  *  Fiv2  lCX  *  CX2  *2  *  CX„  CV  '  V1  c°s  6  cos  * 


V2  [CM  +  CM  e2]  sin 
2  m  Na  Na3  1  V  r 


+  V2 


2,  w  . 


[CNct  +  CNa,  e]  Vsin6  C0S  * 
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.  P  A  pd  V  r„  .  „  2,  w  . 

*  TT-  V  Spa  *  S«3  £  1  V  Sln  * 


+  P  A  Pd  v  r  r  +  c 
4  m  V  luYpa  Sfpa3 


[Spa  *  Spa,  e2)  V  Sin  9  C0S  * 


§-J  v2  CNJ  (‘b  cos  ^  ■  SA  Sin  H  sln  * 
A 


+  y2  Si6  C0S  ^  +  s*n  ^  sdn  6  C0S  ^  =  ® 

A 


C1  "  X0 


r  -  p  A  r 
10  ~  2  m  X2 


C2  “  X0 


C  =  P  A  p 
11  2  m  Sifts 


C3  =  Y0 


c„  =  M-Ed  c 

12  -im  Ypot3 


V  Y0 


^13  =  P  A'  p 
1A  2  m  Xm 


C5  =  Z0 


C14  =  Ss  5A  (2  n) 


C6  *  Z0 


Slfi.  6B  2  n 


c  =  —  ft  c 

7  2  m  X 


C  =  £-A  C 
8  2  m  Net 


C.  =  pd  Cv 
9  4  m  f  Ypa 
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Rewriting  Equation  (43)  with  these  definitions,  one  obtains: 


VXe  +  y2  tC7  +  cio  ^  '  C13  (V  "  V1  C0S  6  C0S  * 
-  V2  [Cg  +  Cn  e2]  (-j-  sin  i|»  -  ^  sin  0  cos  ip) 


+  V  [Cg  +  Cj2  e2]  (y-  sin  ^  +  -■  sin  0  cos  40 


-  V  [(C1c.  cos  <p  -  C  sin  <J0  sin  4*  -  (C1.  cos  <f> 


+  C15  sin  (j))  sin  0  cos  4»]  =  0 


To  perform  differential  corrections  the  3X/3C.  is  required  where  C.  =  the 

•  1  1 

initial  conditions  on  X  ,  Y  Z,X,Y,Z,  and  the  aerodynamic  coefficients 

6  6)  6  6  6  6  • 

(C7  -  C  ).  The  3X/3C.  is  obtained  through  integration  of  3V  /3C..  For 
computational  purposes  3V  /3C.  is  expressed  in  generalized  form  such 

AC 

that  it  will  be  applicable  to  all  C^'s. 

Taking  the  generalized  partial  of  Equation  (44)  with  respect  to  C.,  the 

J 

following  is  obtained: 


+J  V  (2  C?  +  2  C1Q  e2  +  C13  (3  V  -  2  VQ))  cos  0  cos  *  (45) 
-  2  V  (Cg  +  C^,  1  e2)  (^  sin  4>  -  ^  sin  0  cos  40 


+  (Cg  x  C12  e2)  (^  sin  4>  +  ^  sin  0  cos  40 


-  2  V  [(Cjg  cos  <f>  -  C^  sin  <{>)  sin  <p  -  (Cj4  cos  <j> 


+  C,_  sin  $)  sin  0  cos  4>]  f  3C.  =  -  K.  3C./3C. 

15  J  1  111 


L 


where:  K.  =  the  constant  multiplier  of  the  specific  C.  from  Equation  (44), 

1  2  v  w 

for  example,  if  C^.  =  Cg,  then  =  -  V  (^  sin  4»  -  y  sin  0  cos  i{0 


It  should  be  noted  that  v/V  and  w/W  are  trigonometric  functions  of  the 
missile  angles  and  are  assumed  to  be  independent  variables.  0  and  ^  are 
functions  of  0m,  ,  VXg,  Vye,  and  VZe.  However  as  6ffl  and  are  considered 

known  and  the  angular  components  y  and  6g  small,  0  and  ij>  have  also  been  desig¬ 
nated  independent  variables.  It  should  be  pointed  out  that  new  profiles  of 
0  and  are  computed  for  each  iteration  thus  minimizing  this  assumption  pro¬ 
viding  convergence  is  being  achieved. 

3V/8Cj  must  be  expressed  in  terms  of  3VXe/3C^  and  is  done  as  follows: 


2  2  2 

V  s  fv  +V  +  V  1 
v  lvXe  vYe  vXe  J 

Taking  the  generalized  partial  derivative,  one  obtains: 


.  1  ‘2  VXe  SV8Ci  *  2  VYe  8V9Cj  *  2  VZe  9V3Cj1 

av/3Cj  '  2  - - 

[V  2  +  v  2  +  v  2]^ 

L  Xe  Ye  vZe  J 


Vv  3Vv  Vv  3Vv  V-  3V_ 

3V/3C  -  r  +  ^e  .  +  l 

3V/3L.  -  [  v  ac.  V  3C.  V  3C.  1 

J  ill 


Substituting  Equation  (46)  into  (45) ,  the  generalized  partial  derivative 
may  be  written  as  follows: 


3Vv  Vv  3Vv  Vv  3Vv  V7o  3V_ 

Xe  ,  ,  r  Xe  Xe  ,  Ye  Ye,  Ze  le_  ,  _  „ 

3C,  1  1  V  3C.  V  3C .  V  3C.  J  1 

1  3  J  J 


where:  Aj  =  V  (2  C?  +  2  C1Q  e  +  C13  (3  V  -  2  VQ))  cos  0  cos  i|> 


-  2  V  (Cg  +  Cn  e2)  (^  sin  <|>  -  sin  0  cos  *) 

+  (Cg  +  C12  e2)  (~  sin  <J>  +  ^  sin  0  cos  i|») 


-  2  V  l(Cj5  cos  t}>  -  sin  4>)  sin  ^ 

-  (C14  cos  +  C15  sin  f>  sin  0  cos  >|>] 
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1L  =  -K. 
1  J 


Consider  the  remaining  two  equations  of  motion: 

Vye  +  V2  [C?  +  C1Q  e2  +  C13  (V  -  VQ)]  cos  6  sin  ip  (48) 

+  V2  [Cg  +  e2]  (^  cos  ip  +  ^  sin  6  sin  <J>) 

-  V  [Cg  +  C12  e2]  (^  cos  ip  -  sin  0  sin  tf) 

2 

+  V  t(C15  cos  sin  4>)  cos  +  (C^  cos  <j> 

+  C15  sin  sin  0  sin  ij>]  =  0 


VZe  -  V2  [C?  +  C1Q  e2  +  C13  (V  -  Vo)]  sin  0  +  V2  (Cg  +  Cu  e2]  ^  cos  0(49) 

2  v  2 

+  V  [Cg  +  C12  e  ]  Y  cos  0  +  V  (C14  cos  $  +  C15  sin  <j>)  cos  0  -  g  =  0 


Similarly,  as  before,  the  generalized  partial  derivatives  of  Equations 
(48)  and  (49)  may  be  expressed  in  the  following  form: 


3Vv  Vv  3Vv  Vv  3Vv  V_  3V_ 

Ye  ,  ,  r  Xe  Xe  ,  Ye  Ye  ,  Ze  Ze  i 

3?T7"  A2  V  3U7"  V  3C .  V  3C7" 

J  1  1  J 


3V_  VY  3V„  Vv  3Vv  V-  3V_ 

Ze  ,  *  r  Xe  Xe  Ye  Ye  Ze  Ze  i 

3C.  A3  V  3C .  V  3C .  V  3C. 

1  111 


where:  A2  =  V  (2  C?  +  2  C1Q  e*  +  CJ3  (3  V  -  2  VQ))  cos  0  sin  i|> 


+  2  V  (Cg  +  Cn  e2)  (~j  cos  i|>  +  y-  sin  @  sin  <J») 
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2  w  v 

-  (Cg  +  C12  e  )  C^-  cos  i <l>  -  y  sine  sin  ip) 

+  2  V  [(Cjg  cos  <j>  -  sin  41)  cos  ip 
+  (C^4  cos  4>  +  sin  <}>)  sin  0  sin  i/>] 

B2  =  -K.  (of  Vye  equation) 

2 

A3  =  -  V  (2  C7  +  2  C1Q  e  +  C13  (3  V  -  2  Vq))  sin  9 

+  2  V  (Cg  +  Cn  e2)  ^  cos  0 
+  (Cg  +  C12  e2)  ^  cos  0 
+  2  V  (Cj^  cos  +  Cj5  sin  $)  cos0 

B3  =  -Kj  (of  V2e  equation) 

(2)  Derivation  of  Parametric  Differential  Equations  for  Approximate 
Six-DOF  Equations 

Consider  equation  (32) : 

i|)  cos  0  -  M  sin  tp  -  M  ip  cos  0  - 
ma  mq  Y 

"  Mm6^  "  (e  +  8  cos  9/V^ 

-  2  sin  0  (0  +  g  cos  e/V 


Mnpa  sin  0  cos  *  (52) 


^  =  0 
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As  previously  stated,  the  terra,  g  cos  0/V,  is  considered  a  constant  in 
taking  the  partial  derivative  of  Equation  (52) . 


ip  cos  9  -  M  sin  ip  -  M  ip  cos  0  -  M  sin  0  cos  ip  (531 

ma  raq  npa  v 


■  MmSip  -  (e  +  K  )  p  r  -  2  sin  6  (0  +  K_ )  *  =  0 


where:  M  =  9JL1  [c  +  C  e2  +  C  c4  +  C  (VREF  -  V)  ] 

ma  I„  ma  ma3  ma5  ma 


M  =  3_Ai_  [C  +  C  e2  -  C  2  1  /rad21 
mq  2  V  I  mc*  mq2  Not  y 

y 


M  =  §-4-4 — E-  [C  +  C  e2  +  C  e4  +  CKf  2  I  /md2] 
npa  2  V  I  nna  nt>a<»  m>a<;  Na  x'  J 


npa  npa 3 


npa  5 


Mm««  *  V-4  [CmS  (SB  sin  *  -  SA  cos 


Mm«*  *  [Cm6  (SB  cos  *  *  SB  sin  *' 


2  2  —  2  2  2 
e  =  sin  a  =  sin  ip  +  sin  0  cos  ip 


K  -  g  cos  0/V 

o 

To  perform  differential  corrections,  the  3ip/3C.j  is  required  where  Cj  =  the 
initial  conditions  on  0,  <p,  0,  ip,  and  the  aerodynamic  moment  coefficients. 

The  3ip/3Cj  is  obtained  through  integration  of  3ip/3Cj .  For  computational  pur¬ 
poses  3ip/3Cj  is  expressed  in  generalized  form  such  that  it  will  be  applicable 
to  all  C. *s.  The  angular  motion  reduction  follows  the  completion  of  the  roll 
and  translational  motion  reductions  and  thus  the  quantities  V,  q'  and  p  are 
considered  known  and  can  be  designated  independent  variables. 

Taking  the  generalized  partial  of  Equation  (53)  with  respect  to  C.,  the 
following  is  obtained: 
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(54) 
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(56) 


g-  cos  9  -  *  sin  9  f§-  -  um  cos  ♦  ~  -  M  cos  9  g- 
3  3  3  3 


+  M  ip  sin  6  - M  cos  0  cos  ip  +  M  sin  0  sin  ip  .^5— 

mq  3C.  npa  r  3C.  npa  y  3C- 

3  3  3 


-  P  y-  -  2  sin  0  (6  +  K  )  ||-  -  2  sin  0  ip  -  2  cos  0  (6 

+  v  ♦  M: 

*  3 


3M 


3M 


3M 


-  [sin  ip  — 52  +  ip  cos  0  —23.  +  sin  0  cos  ip  — 223.  ]  [cos2  0  sin  2 ip 


3e 


3e 


3e 


dip 

3C. 


'•3 


7 

+  cos  ^  sin  2  0  ]  = 

j 


-  F.  3C./3C. 
3  3  3 


Equation  (56)  may  be  written  in  the  following  generalized  form: 


dip  .  .  30  „  30  „  3ip  .  3i|i 

sc:  cos  6  +  Ai  jzr +  Bi  3^7  +  ci  dir +  Di  z: 

3  3  3  3  3L 


=  Glj 


3 


where:  A,  =  -  ip  sin  0  +  M  ip  sin  0  -  M  cos  0  cos  ip 

mq  npa  y 


' 1 


3M 


-  2  cos  0  (0  +  K  )  ip  -  [sin  ip  — *2.  +  ip  COs  0  — j 


3M 


m 


3e 


3e 


3M 


+  sin  0  cos  ip  ■  —222.]  [cos2  ip  sin  2  0] 

3e 


Bj  =  -  p  2  -  2  sin  0  ip 

y 
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3M 


C,  =  -  M  cos  ip  +  M  sin  0  sin  \p  -  [sin  ip  — ^ 
1  mo  npa  Y  f  2 

o£ 


ma 


3M  3M 

+  tfi  cos  0  — +  sin  0  cos  tp  — 2E a ]  sin  2  ip  cosZ  0 
3e  3e 


D,  =  -  M  cos  0-2  sin  0  (0  +  K  ) 
1  mq  v  zJ 


G, .  =  -  F. 
ll  1 


From  the  definitions  of  M  ,  M  ,  M  ,  the  following  partial  derivatives 

ma  mq  npa  °  r 

are  used  in  Equation  (56) ; 


3M 


3e 


=  sAA  [C  +  2  C  ] 
2  I  ma3  ma5 


3M  2 

. .  HfL  =  qAd  [c  ,] 
8e2  rvr-  “>2 


3M  -  .  ,2 

=  aJLl-z  rc  +  2  c  ] 

2  2  V  I  1  npa*  npa5J 


3e 


Consider  the  remaining  equation  of  motion  used: 


0  -  M  sin  0  cos  ip  -  M  (0  +  K)+M  sin  i[»  +  M  ,Q 
ma  mq  g  npa  m60 


(57) 


*2 


+  ip  cos  0  p  j—  +  ip  cos  0  sin  0  = 

y 


_  0 


Similarly,  as  before,  Equation  (57)  may  be  expressed  in  the  following 
generalized  form: 


90  &  36  *  r  96  *  r  3#  +  n  3*  -  r 

3C.  A2  3C .  B2  3C .  C2  3C.  °2  3C.  "  G2j 

1113  3 
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where: 


•  x  *2 

A0  -  -  P  ^  sin  0  - - M  cos  1 p  cos  G  +  i|>  cos  2  6 

2  I  mot 


3M  .  3M  3M 

-  [  cos  ip  sin  ©  — 2J2  +  0  — 29.  _  sin  ip  — 9E9.  ]  cos^  ^  sin  2  0 

3e  3e  3c 


B0  =  -  M 

2  mq 


CL  =  M  sin  ip  sin  0  +  M  cos  ip 
2  ma  T  npa 


•  an  an  an  n 

[  e  ..29.  +  cos  ip  sin  0  — 29.  _  sin  ip  — -229.  ]  cos  e  sin  2  \ p 


D2  =  p  cos  0  +  sin  2  0  ip 

Y 


G2j  =  -  Fj 


(3)  Derivation  of  Partial  Derivations  for  Approximate  Roll  Equation 


Consider  the  equation: 


q  A  d2 


♦  '  Kt-  tClP  +  Clpv  <Vo  -  <♦  -  *  Sin  9>  -  H^C16  6 

XV  X 


-  tan  0  (tp  cos  0  -  ip  0  sin  0)  + 


cos  9 


The  experimental  data  is  provided  as  the  angle  <p. 


The  initial  conditions  are  <j>Q  and  <Pq. 


M  ^  2  • 

jL  1  _  9  A  d  3  4>  rr  +  r  fV  -  Vll 

3  C.  2  V  I  3  C.  lClp  Clp,,  tVo.  V)J 

2  x  j  r  rV 


-  F.  3C./3C. 

j  r  i 


41 


Use  is  made  of  the  equality  p  =  $  -  ip  sin  0.  0,  ifr,  i|/,  V  and  q  are 

considered  known  and  therefore  are  designated  independent  variables.  This 
reduces  the  complexity  of  this  problem  as  most  higher  order  terms  (0,  t|>,  <j>) 
now  vanish.  The  spin  (p)  was  computed  during  the  integration  of  the  equations 
of  motion,  and  as  C4>  -  sin  0)  always  appears  as  a  paif,  p  is  substituted. 
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SECTION  V 

TECHNIQUE  APPLICATION 


1.  COMPUTER  PROGRAMS 

The  equations  derived  in  the  previous  section  are  applied  through  five 
computer  programs  interconnected  through  the  use  of  overlays. 

The  five  programs  consist  of: 

1.  PRENUT  -  A  modified  linear  theory  program  used  for  initialization. 
This  program  provides  sufficient  initial  guesses  and  program  con¬ 
trols  to  properly  execute  the  numerical  integration  data  reduction 
process . 

2.  EARTH  HEEVE  -  This  program  analyzes  the  translational  and  roll 
motion  of  the  projectile,  determining  the  earth-fixed  velocity 
vectors  and  roll  rate.  It  transforms  the  measured  fixed  plane 
angles  to  missile  angles  and  determines  the  axial  force,  normal 
force  and  roll  coefficients. 

3.  MODYAW  -  This  program  makes  a  final  complete  modified  linear 
theory  analysis  of  the  rolling,  translational,  and  angular  motion 
of  the  projectile  from  the  experimental  data. 

4.  TRIANGLES  -  This  numerical  integration  program  is  utilized  to 
analyze  the  angular  motion  of  one,  two,  or  three  experiments 
simultaneously.  Several  options  are  available  to  the  user  and 
can  be  utilized  for  obtaining  non-linear  aerodynamic  coefficients. 

5.  TRIHEEVE  -  The  translational  motion  of  up  to  three  experiments  is 
fitted  by  this  program.  Options  are  available  to  assist  the  user 
in  obtaining  non-linear  aerodynamic  coefficients. 
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Each  program's  inputs  and  outputs  except  MODYAW  are  discussed  in  the 
following  section.  A  diagram  showing  pictorially  the  flow  of  a  typical 
reduction  of  three  sets  of  data  which  are  to  be  reduced  simultaneously  by 
TRIHEEVE  and  TRIANGLES  is  given  as  Figure  3. 


DIRECTION  COSINES 
EXPERIMENTAL  xfi  y#  z# 
EXPERIMENTAL  ROLL  ANGLE 


After  three  designated  experiments  are  reduced  individually, 
a  multiple  TRIHEEVE  reduction  is  made  followed  by  a 
multiple  TRIANGLES  reduction. 


Figure  3.  Data  Reduction  System  Flow  Diagram 
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a.  Program  -  PRENUT 


Purpose  -  Perform  a  linear  theory  data  reduction  and  provide  estimates  of 
initial  conditions  and  aerodynamic  coefficients  for  the  EARTH 
HEEVE  program  and  the  TRIANGLES  program. 


Inputs 


-  Earth  fixed  direction  m  ,  n  ,  p 

e  e  e 


cosines 

Earth  fixed 
coordinates 

Roll  angle 


X  ,  Y  ,  Z 
e’  e  e 


>  versus  t 


Projectile  physical  m,  I  ,  I  ,  CG,  d,  a 

x  y 

properties 

Estimates  for  C„ 

Ncs 


Outputs  -  Experimental  data 
Fitted  values  of 


'J,»  4* >  Xe,  Yg,  Ze,  t 

0  ip  versus  t 
m,  m 


Initial  conditions 


X  ,  X  ,  Y  ,  Y  ,  Z  ,  Z  versus  t 

eo  eo  eo  eo  eo  eo 


«  • 


\>  %■  *0‘  +o>  Po 


Aerodynamic 

coefficients 


Cvn,  C  ,  C  ,  C  ,  Cl, 
XCr  ma'  mq'  npot'  Tla 


Comments:  The  roll  angle,  <f>,  is  measured  from  0°  to  360°  degrees  and  must 
be  correctly  (exactly)  unwound  to  be  a  continuous  accumulative 
function  of  time  (0°  to  «>) .  This  is  accomplished  in  a  subroutine 
contained  in  PRENUT. 
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Equation  of  Motion  -  PRENUT 


K  (*1  +  iw1)Xe  +  (X2  +  iw0)X 


2'  e  *  K. 


amplitude  nutation  vector 
K2  amplitude  precession  vector 

Aj  nutation  damping  factor 

>>2  precession  damping  factor 

nutation  frequency  (<{>  +  <p  X  ) 

x  x  6 

_  •  •• 

precession  frequency  (<p  +  <f>  X  ) 

yaw  of  repose  vector 
Xe  distance  traveled 

nutation  vector  orientation 
$2  precession  vector  orientation 


< 


1 


I 

1 


i 


1 


4 


b.  Program  -  EARTH  HE EVE 


Purpose  -  Perform  a  numerical  integration  reduction  of  the  earth- fixed 
(translational)  equations  of  motion  (equations  40,  41  and  42) 


determining  best  fit  to  X  ,  Y  ,  Z  . 

0  6  6 


Inputs 


-  From  PRENUT 
Experimental  data 
Fitted  data 
Physical  properties 
Estimates  for 


0,  '!»,  4>>  Xg,  Yg,  Ze,  versus  t 


0  ,  versus  t 
m  m 


m,  I  ,  I  ,  CG,  d,  i 
*  x  y 


cvn,  cxt  ,  cv_,  cv  ,  c  ,  c  ,  c  ,  c, 

XO*  Net*  X2*  Ypa  ma’  npa  mq  lp 


Initial  conditions 


X.  »  XP  >  Ye  *  Ye  »  Ze  »  Ze  *  VerSUS  1 
e0  0  0  0  0  0 


•  • 


6m0*  *m0»  V  V  V  p0 


Outputs  -  Experimental  data 
Computed  data 


©,  ♦  »  *e»  ye,  ze  versus  t 


0  ,  ^  versus  t 
nr  m 


Computed  initial 
conditions 


X  ,  X  ,  Y  ,  Y  ,  Za  ,  Za 

eo  eo  eo  eo  eo  e0 


Computed  parameters 


V,  q,  p,  yg  versus  t  (from  best  fit) 


Computed  coefficients  Cx0,  C^,  Cypo>  Cxv,  Clp, 


Comments:  The  final  numerical  integration  roll  fit  is  performed  in  this 


program  with  final  values  of  C^  and  C^  determined.  Exact 


transformation  from  mfeasured  fixed  plane  angles  to  missile  angles 


are  performed  using  the  theoretical  X  ,  Y  ,  and  Z  velocity  profiles. 

6  6  6 
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3[>  3 1  >  3 1>  3|>  3|>  3 1  >  3 1  >  3|>  3|>  3|>  3|>  3|> 


Equations  of  Motion  -  EARTH  HEEVE 


[C^  +  C^2  e  3  cos  0  cos 

[CM  +  CM  e2]  sin  t|» 

Na  Na3  V 

[CN«  *  V  ^  ?  sin  6  cos 

E-i  tcv  *  Cv  e2]  £  sin 
2  V  Ypa  Ypa3  V 

?-  y.~  [Cv  +  Cv  e  ]  tt  sin 
2  V  Ypa  Ypc»3  V 


[Cx  +  CX2  e  ]  cos  0  sin  \p 

[CM  +  CM  e2]  cos  4> 

Na  Na3  V 

2  w 

[C.,  +  CM  e  ]  rr  sin  0  sin 
Na  Na3  V 

2-4  [Cv  +  Cv  e2]  tt  cos 
2  V  Ypa  Ypa3  V 

2-4  [Cy  *  Cv  e2]  4  sin 
2  V  Ypa  Ypa3  V 


[Cx  +  Cx2  e  ]  sin  0 


(CNa  +  CNc.3  *2]  V  005  6 


qAE_i[r  +  Cv  e2]  J  cos 
n  u  2  V  Ypa  Ypas  V 


P .  Vr^  'V *  "V 

x 


1PJ 


4> 

0  COS 


* 

<1» 

0  sin  <J> 


0  -  g 
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c.  Program  -  TRIANGLES 


Purpose  -  To  provide  a  best  fit  to  the  angular  motion  experimental  trans¬ 
formed  data  (0^,  iji  )'  using  the  equations  of  motion  (equations  30 
and  31) .  To  provide  a  multiple  fit  of  more  than  one  experimental 
case  using  common  coefficients  for  all  cases. 


Inputs  -  From  EARTH  HEEVE 


Experimental  data 

e»  xe>  Ye>  ze»  versus 

Transformed  data 

0  ,  t  ,  versus  t 
m  m 

Computed  data 

V,  9*  P>  versus  t 

C 

Computed  coefficients 

CX0'  CNa’  Spa"  CXV>  CX2 

Estimated  coefficients 

C  ,  C  ,  C 
ma  mq  npa 

Estimates  for 

•  • 

0  ,  <p  ,  0,  ij> 

ra0  0 

Physical  properties 

m,  I  ,  I  ,  CG,  d,  a 
x  y 

Initial  conditions 

Xe  ’  Xe  *  Ye  ’  Ye  »  Ze  ’  Ze 
e0  e0  0  e0  e0  e0 

Outputs  -  Experimental  data 
Best  fit 
Computed  data 

Computed  coefficients 

Final  computed 
coefficients 
Initial  conditions 


0,  tf>,  <J>,  Xe,  Ye,  Ze,  versus  t 

0  ,  versus  t 
nr  rm 

Pr  versus  t,  I  /I 
x  y 

CX0’  Sla’  ^Ypct’  CXV’  CX2 


C  ,  C  C  plus  non-linear  terms 
ma  mq,  npa  r 


0 
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Comments:  The  final  aerodynamic  moment  coefficients  Cma>  C  ,  and 

(functions  of  angle  of  attack  and  Mach  number)  are  derived  from 
this  program  both  for  single  and  multiple  fits. 


( 


Equations  of  Motions  -  TRIANGLES 


u=-q  —  Cy  +  g  sin  6  -  0  w  +  t  cos  6  v 


V  =  -  o'  -  CL  —  +  q  —  EjL  Cv  ~  -  <P  cos  e  [u  +  w  tan  6]  (2) 

v  q  m  Tla  V  4m2V  Ypa  V 

w  -  -  a  -  C.  -  -  Q  -  C„  ^  -  g  cos  6  +  <J»  sin  0  v  +  6  u  (3) 
w  q  m  la  V  "  m  2  V  Ypa  Vs 

•  •  •  •  •  • 

a  =  a  LAe£  r  +  tan  6  [cos  0  <p  -  ip  0  sin  6]  +  6  ip/cos  6  (4) 

4I  2  V  lp 


p  =  <j>  -  i|»  sin  0 


S.flAltC  +  C  e2  +  C  e4  +  C  e6*C  CV0-WjT  C5) 

0  q  I  1  ma  ma3  mas  ma7  raav  0  v 


+  r  LA —  [c  +  c  e2]  0 

q  2  V  I  1  mq  rnq9 

V 


—  A  a  p  rr  r  2  c  e^i  X. 
+  q  2  V  I _  Cnpa  Cnpa3  npa5  V 


-qr^Sns  (6B  Sin*-6A cos  *)] 

y 

1  *2 

-  tji  cos  0  p  j— ■  -  ^  cos  0  sin  0 

y 


■[ ' -  <  f4  *  C»a,  ^  +  C».5  **  *  C-7  +  C-v  (V«>f 


K  A2  2  • 

+  qLA—  [C  +  C  G  3  <f>  COS  0 
q  2  V  I  mq  mq2 
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+  q  2  V  1^  CCnpa  +  CnP“3  £  +  CnP«5  £  ^  v 


+  ?  tCm6  (6B  C0S  ♦  *  6A  Sin  « 1 
y 


•  •  •  Ix1 

+  2  8  f  sin  0  +  ®  P  — jcos  6 


d.  Program  -  TRIHEEVE 


Purpose  -  To  provide  a  best  fit  of  the  equations  of  motion  to  the  measured 
X,  Y,  and  Z  data.  To  provide  a  multiple  fit  of  more  than  one 

experimental  case  using  unique  coefficients  for  all  cases. 

Inputs  -  From  TRIANGLES 

Experimental 

Q,  Xe,  Ye>  Ze,  versus  t 

Best  fit 

9  \b  versus  t 
m,  rm 

Computed  coefficient 

CX0*  W  CYpa>  CX2*  CXV 

(estimates) 

,  ♦  • 

Computed  initial 
conditions 

x  ,  X  ,  Y  ,  Y  ,  z  Z 
eo  e0  e0  eo  eo  eo 

Computed 

-Pr  vs  t 

Physical  properties 

m,  Ix>  Iy,  CG,  d,  Z 

Outputs  -  Final  best  fits  of 

X  .  Y  ,  Z  versus  t 
e’  e*  e 

Computed  data 

V,  q,  p,  Ye  versus  t 

Final  computed 

r  C  C  as  function  of  Mach  number 

WX’-  Na’  Ypa 

coefficients 

and  angle  of  attack 

Transformed  0m  and  <|<m 

based  on  new  velocity  profile 

Comments:  TRIHEEVE  computes  the 

final  best  value  of  the  axial  force, 

normal  force,  and  Magnus  Force  coefficients  both  for  single 
and  multiple  fits. 
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Equations  of  Motion  -  TRIHEEVE 


vXe  =  -  “‘[Cx  +  CX2  sin2  a  +  CXV  (V  '  VREF^  cos  6  cos  * 

*  [CNo  *  V  s1"2  lV  sln  *  '  7  sln  6  cos  I11 

-  3iVvi  V  +  Vs  si"2  sin  *  *  7  cos  *  sin  61 

Q  A 

+  -! —  [-C...  (sin  <i  sin  \p  +  cos  <{>  sin  6  cos  </>) 
m  no, 

A 


+  C...  (cos  <J)  sin  ip  -  sin  <j»  sin  6  cos  i|»)] 

N6B 

Vv  =  .  1A  [c  +  C  sin2  a  +  CY..  (V  -  VREF)]  cos  6  sin  ip 
Ye  m  X  \2  av 

—  2 

-  5LA  [C  +  CM  sin  a]  cos  ip  +  ^  sin  6  sin  U»] 

m  Na  Na3  V  V 

+  ^  tcv  +  Cv  sin2  o]  [rr  cos  ip  -  sin  ^  sin  0] 

m  2  V  Ypa  Ypa3  V  V 

-  S_A  [_  r  (sin  <p  cos  ^  -  cos  <f>  sin  0  sin 

m  N6.  v 
A 

+  C...  (cos  4>  cos  +  sin  «J>  sin  0  sin  ij/)] 

N0_ 


*  V  [CX  *  CX2  s1"2  “  *  CXV  <V  -  VREF)1  sin  9  -  *r  [CNa 
2  w 

+  C.,  sin  ]  tt  cos  0 
N(*3  V 


.  a.  A-iLi!.  i 

c  + 

c 

hTTT  1 

1  Yptx 

Yjm., 

+  CMX.  sin 
N6b 

<j )]  cos 

e  -  g 

.2  ,  v 

Sin  it]  y  cos  0 


q  A 


[C  cos  * 
A 
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2.  NUMERICAL  ROUTINES 


All  numerical  integration  performed  in  the  previously  described  programs 
utilize  the  fourth  order  Runge-Kutta  method.  It  is  acknowledged  that  there 
are  extensive  numerical  integration  techniques  in  existence,  many  of  which 
are  claimed  to  be  much  more  efficient  than  the  Runge-Kutta  method.  The  two 
primary  reasons  for  implementing  this  technique  are  that  (1)  Runge  Kutta 
methods  are  self-starting,  the  interval  between  steps  may  be  changed  at  will, 
and  in  general,  they  are  particularly  straightforward  to  apply  on  a  digital 
computer,  and  (2)  they  are  comparable  in  accuracy  and  often  more  accurate 
than  corresponding  order  predictor-corrector  methods  except  that  the  integra¬ 
tion  interval  is  a  more  significant  parameter  in  terms  of  solution  accuracy. 
Experience  has  demonstrated  this  technique  to  be  reliable,  stable,  and  of 
sufficient  efficiency  from  a  storage  versus  execution  time  viewpoint  given 
the  form  of  the  differential  equations  solved  for  as  previously  described 
(equations  of  motion  and  partial  differential  correction  equations) .  During 
development  of  these  programs,  test  implementation  of  a  modified  Adams-Moulton 
predictor- corrector  routine  was  made  for  purposes  of  comparison.  The  results 
indicated  that  both  integration  techniques  achieved  the  same  degree  of  accur¬ 
acy  and  execution  time.  Comparable  execution  times  resulted  from  the  fact 
that  the  predictor-corrector  method  had  difficulty  in  starting  and  time  step 
adjustment. 

The  matrix  inversion  routine  implemented  in  all  the  programs  utilizes  a 
Gaussian  elimination  method.  This  is  a  direct  method  (versus  iterative)  and 
eliminates  truncation  error.  Over  a  period  of  years  during  which  many  sets 
of  test  data  were  reduced,  no  problems  were  encountered  with  either  the 
numerical  integration  or  the  matrix  inversion  methods. 

Additional  details  on  these  two  numerical  routines  may  be  found  in 
references  26  and  27. 
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SECTION  VI 

CONCLUSIONS  AND  FUTURE  ACTION 


An  extensive  set  of  data  reduction  computer  programs  for  analyzing  aero- 
ballistic  range  data  has  been  developed  and  sequenced  to  produce  a  sophisti¬ 
cated  system.  These  programs  are  primarily  geared  to  the  reduction  and 
analysis  of  spin- stabilized  or  spinning  statically  stable  projectiles  with 
small  configuration  asymmetries. 

The  need  exists  for  a  similar  set  of  programs  tailored  toward  the  analy¬ 
sis  of  rolling  or  non-rolling  projectiles  with  configurational  and  mass 
asymmetries.  In  this  manner,  the  ballistic  range  could  become  a  useful  tool 
in  the  testing  of  winged  re-entry  bodies  and  non-symmetric  nose  cones. 


DOWN  RANGE 


SIDE  VIEW 

POSITIVE  DIRECTIONS  SHOWN 

Figure  1-1.  Coordinate  System  Fixed  Plane  (Positions) 
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Angular  Identities 
Missile  Angles 


sin  =  V  sin2  +  cos2  ^  sin2  ^ 

si"  8b  ‘  sin  \/cos  *m 

sin  em  =  sin 

sin  ♦„  -  '  sin  6,  cos  xn 

tan  6m  =  -  tan  *m/cos  0m 

Missile  Angles  -  Earth  Angles 

0  =  0  -  sin"1  [sin  y  /cos  ] 

m  e  m 

ijj  =  sin"1  [sin  ^/cos  Yel  +  6e 

Angular  Identities 
Earth  Angles 


tan  6  =  tan  x  cos  <J> 


Missile  Velocities  -  Missile  Angles 

v  =  -  V  sin 

w  =  V  cos  <|>  sin  9 
m  m 

u  =  V  cos  ^  cos  em 
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Missile  Angle  -  Velocity  Relationship 


ANGLE 

SINE 

' 

COSINE 

TANGENT 

■ 

-y  2  2 

v  U  +  V 

w 

w 

,/  2  2  2 

J  2  2  2 

J  2  2 

V  U  +  V  +  w 

V  U  +  V  +  w 

V  U  +  V 

V  2  2 

K 

-V 

V  u  +  w 

-V 

in 

J  2  2  2 

J  2  2  2 

J  2  2 

yu  +  v  +  w 

Vu  +  v  +  w 

VU  +  w 

w 

u 

w 

e 

m 

/  2  2 

Vu  +  w 

J  2  A  2 

vu  +  w 

u 

a 

m 

w 

u 

w 

/  2  2 

/  2  2 

u 

-yu  +  w 

•yu  +  w 

V 

u 

V 

— 

l/u2  *  V2 

Vu2  +  V2 

u 

,/  2  2 

Vv  +  w 

u 

a 

m 

-i/  2  2  2 

u 

V  u  +  v  +  w 

■GSM 

Trajectory  Angle  -  Velocity  Relationships 


ANGLE 

SINE 

COSINE 

TANGENT 

VZe 

VZe 

Ye 

Vv  2*v  2 

v  vXe  vYe 

6 

^MB 

VYe 

e 

Bms 

VXe 

Transformations 
Missile  to  Earth 

VXg  =  u  cos  e  cos  tj>  -  v  sin  ip  +  w  sin  0  cos  \p 

Vye  =  u  cos  6  sin  ij/  +  v  cos  ip  +  w  sin  0  sin  ip 

V-  =  -  u  sin  0  +  w  cos  0 

Ze 


Direction  Cosines  to  Fixed  Plane  Angles 

6  =  -  sin-1  (n) 

¥  =  sin_1(ine/-\/me2  +  pe2) 

where:  n  ,  m  ,  p  are  the  direction  cosines  of  the  missiles  X  axis  relative 
e  e  re 

to  the  range  coordinate  system  (fixed  plane) . 
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